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Abstract 

Dispersive constraints on the shape of the form factors which describe the exclusive 
decays B D^*H V are derived by fully exploiting spin symmetry in the ground- 
state doublet of heavy-light mesons. The analysis includes all twenty B^*^ D^*^ 
semileptonic form factors. Heavy-quark symmetry, with both short-distance and 
1/mQ corrections included, is used to provide relations between the form factors 
near zero recoil. Simple one-parameter functions are derived, which describe the 
form factors in the semileptonic region with an accuracy of better than 2%. The 
implications of our results for the determination of \ Vcb\ are discussed. 
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1 Introduction 



The differential rates for tlie exclusive semileptonic decays B D'^*Hu are given by [|l| 



1 + 



{niB — mD*Y — l{w + 1)^ 

l^(^) 



dw 487r^ 

Aw m"^ — 2w msmj:,* + m'jj. 



w + 1 {niB — mD*Y 

"^'''7/''' = ^|||^(™B + m.)',4(»^-lFmWP. (1) 

where J^{w) and Vi{w) are hadronic form factors, and w = vb ■ v^m is the product of 
the velocities of the initial and final mesons. In the heavy-quark limit, J-'{'w) and Vi{w) 
coincide with the Isgur-Wise function ^{w), which describes the long-distance physics 
associated with the light degrees of freedom in the heavy mesons This function is 
normalized to unity at zero recoil, corresponding to w = 1. Corrections to this limit can 
be calculated using the heavy-quark effective theory [Q, and are suppressed by powers of 
asimq) or Aqcd/^^q? where we use ttiq generically for or rric. Detailed calculations 
of these corrections lead to = 0.91 ± 0.03 [|, and Vi{l) = 0.98 ± 0.07 g], so 

that an accurate determination of \Vcb\ can be obtained by extrapolating the differential 
decay rates to w = 1. To reduce the uncertainties associated with this extrapolation, 
constraints on the shape of the form factors are highly desirable. A suitable framework 
to derive such constraints is provided by a dispersion technique proposed some time ago 
H^EMl' ^^"^ been applied to heavy-meson form factors more recently H-ii. This 



method is based on first principles: the analyticity properties of two-point functions of 
local currents and the positivity of the corresponding spectral functions. Analyticity 
relates integrals of the hadronic spectral functions to the behavior of QCD two-point 
functions in the deep Euclidean region via dispersion relations. Positivity guarantees 
that the contributions of the states of interest to the spectral functions are bounded 
above. The constraints on the relevant form factors, given their analyticity properties, 
then follow from these bounds. 



Following Refs. |T^, ^ , we generalize this method to fully exploit spin symmetry in 
the doublets of the ground-state B^*^ and D^*^ mesons. By taking into account all contri- 
butions related by heavy-quark symmetry, this generalization increases the constraining 
power of the method considerably. We investigate all spin-parity channels (J^ = 0"^, 0~, 
l" and 1+) relevant for B^*^ D^*^ transitions. This gives us four inequalities and thus 
four allowed domains for the parameters which describe the form factors of interest. The 
optimal constraints are obtained by taking the intersection of these domains. Although 
these different channels were also investigated recently by the authors of Ref. |2^, we 
bring a number of improvements to their analysis and exploit the dispersive techniques 
with a somewhat different emphasis. First, we include the contributions of the B^. poles 
to the relevant polarization functions. Since these contributions are positive, they reduce 
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the domain available to the contributions of the B^*^D^*^ states and therefore increase the 
constraints on the corresponding form factors. As a second improvement, we include the 
short- distance and l/rriQ corrections to the heavy-quark limit in the relations between 
form factors, thereby eliminating the leading uncertainties in these relations. We further 
discuss the role of remaining uncertainties and show how their inclusion is necessary 
to avoid overinterpreting the results of the dispersive method. In addition, instead of 
requiring rather complicated kinematical functions and the knowledge of the singularity 
structure of individual form factors, we introduce single-parameter descriptions that are 
simple, low-order power-series expansions in kinematical variables, and accurate to 2% 
over the full semileptonic domain. Thus, these parametrizations will be very useful for 
reducing the uncertainty in the extraction of \Vcb\ from experimental data. 

The paper is organized as follows. In the next section, we present the derivation of 
the unitarity inequalities, which are the basis of the dispersive bounds on form factors 
discussed in Section ^ In Sections | and |^, we derive model-independent bounds on the 
heavy-meson form factors using expansions in different kinematical variables. We present 
a new method for including corrections to the heavy-quark limit, which fully takes into 
account the uncertainties in their calculation. We find that these uncertainties have 
a significant effect on the bounds and, ultimately, limit the accuracy of the dispersive 
approach. We also derive simple, one-parameter descriptions of the form factors in the 
semileptonic region. The phenomenological consequences of our results are discussed in 
Section ^ where the reader not interested in the technical details of our work can find a 
summary of our parametrizations for the B D^*^^ form factors. Our conclusions are 
given in Section ^ The paper also comprises an Appendix, where we give the definitions 
of form factors and details of our calculations. 

2 Dispersion relations and unitarity inequalities 

We start with the vacuum polarization tensor 

W{q) =i yd^xe*'?-^(0|T{V^(x),V^'^(0)}|0) 

= (gV-^7^V)ni.(g2)+gVno+(g'), (2) 

and a similar expression for the correlator of two axial currents with invariant functions 
ni+ and IIq-, where the subscripts indicate the spin-parity quantum numbers of the 
intermediate states contributing in the various channels. The invariant functions satisfy 
subtracted dispersion relations, which we write as {Q^ = —q^) 
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In the Euclidean region, where > 0, these functions can be calculated in QCD using 
the operator product expansion. On the other hand, the spectral functions Imno±(t) 
and Imni±(t) are given by unitarity relations. In the case of the vector correlator, we 
have 

{qf^q" - g^"'q'^) ImHi- (g^ + ie) + g'^g'^ Imno+ (g^ + ie) 

= lf^dpr (2vr)^5(^)(pr - g) ( |r'^(0)|r) (r|Vt^(0)| ) , (4) 
r 

and a similar relation holds for the axial correlator. The spectral functions, which are 
sums of positive terms, are bounded below by the contributions of the two-particle 
states \BD), \BD*), \B*D), and \B*D*), which are related through crossing symmetry 
to the matrix elements relevant for semileptonic decays. Combining this lower bound 
with the QCD result for the correlators, we derive constraints on the form factors in 
the semileptonic region. These constraints are further improved by using heavy-quark 
symmetry to relate the matrix elements associated with these different contributions. 
Note that, when extending the number of hadron states in the unitarity sum, there is, 
in principle, a danger of double counting if particles which decay via strong interactions 
are included. In our case, D* accounts partially for the two-particle state Dtt, which 
therefore must not be included separately. {B* is stable with respect to strong decays, 
since uib* < + m.„.) 

To further increase the constraining power of the dispersive method, we take into 
account the single-particle contributions of the first two vector B* and pseudoscalar B^ 
mesons.^ These contributions are determined in terms of the masses and leptonic decay 
constants of these states, whose values can be calculated rather accurately using quark 
models. In Table |1|, we collect the masses and decay constants of all Be mesons below 
the threshold for B*D* production, as predicted by the model of Ref. Note that the 



scalar and axial states, as well as some of the vector and pseudoscalar states, have non- 
vanishing orbital angular momentum and therefore vanishing decay constants. Although 
this conclusion is strictly only valid in the context of quark models, we expect orbitally 
excited states to have much smaller decay constants than S-wave states, and take the 
conservative approach of neglecting their contributions to the various spectral functions. 

Inserting these one- and two-particle states into the unitarity sums, and combining 
the resulting bounds with the dispersion relations of (|^), we find 



dw{w -1)^ [^^] {i + Sj2)—. < — Xo+ 



^We do not include more massive resonances, because they lie above the BD threshold and may lead 
to double counting. 
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Table 1: Masses (in GeV) and decay constants (in MeV) of the lowest-lying Be 
states 
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+ E /d«;(^^-l)^/^(^) (1 + 5,7 

i=5,6,7 



2 _ w+l Y 

i 2 y 



< mB*mD* Xi+ (0) , (5) 

^/ 

where we have set = for simphcity. In principle, the inequalities could be strength- 
ened by using negative values of Q^, which however must be sufficiently far away from 
threshold in order for the operator product expansion to be well defined. However, in 
Ref. |jl9[ it was found that this would improve the bounds only slightly. Therefore, we 
shall not consider this possibility any further. 

In terms of the kinematical variable w entering the unitarity inequalities, the mo- 
mentum transfer t = appearing in and (|^) is given by 

t = = m|(.) + m^(.) - 2m^(.)m£,(,)ti? . (6) 
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We work with w instead of t because it is the natural variable for implementing heavy- 
quark symmetry, and because the various B thresholds, which occur at different 
values of t, all occur at the same value w = —1, enabling a unified treatment of all form 
factors. The form factors Vi, Si, Ai and Pi, which appear in @, are linear combinations 
of the traditonal heavy-quark basis form factors chosen (i) to reduce to the Isgur-Wise 
function in the heavy-quark limit, (ii) to have definite spin-parity quantum numbers 
0+, 1+ and 0~, respectively), and (iii) to reduce the unitarity relations to sums of 
squares. The definitions of these form factors are given in the Appendix. The factor of rij 
appears in (|^) because SU(n/) light-flavour multiplets of heavy- meson states contribute 
with the same weight to the unitarity sums. To be conservative, we take Uf = 2.5 to 
account for the breaking of SU(3) flavour symmetry due to the strange-quark mass. For 
each B^*^ D^*^ transition form factor, we define a quantity 

where it is understood that the appropriate masses are used. Moreover, we use the 
notation 



Xi-(O) =Xi-(0)- E 



t2 KiB*) ' 

Xo-(0)=Xo-(0)- E S§-y (8) 

n=l,2 J-^^n\.Bc) 

and thus include the contributions of the first two vector 5* and scalar mesons to 
the polarization functions. 

For the vector correlator, the leading terms in the operator product expansion of the 
vacuum polarization functions are given by 

2n 3 /• / xml + {I - x)ml - m^rric 



Att'^ j xml + (1 — x)m'^ + x{l — x)Q'^ ' 

2x 3 ] 2 3a;m^ + 3(1 - x)m^ + m;,mc + 2x(l - a;)(52 

Ki-yQ ) = — ^ dxx (I — x) 9 , 9 

87r2i ^ ' [xml + {1 - x)ml + x{l - x)Q^f ^' 

where rrib and are the heavy-quark masses. The functions Xo- ^^id Xi+ for 
axial correlator are obtained from these expressions by the replacement rric — > —rric. 
Evaluating these results at = 0, and using z = mc/rrib = 0.29 ± 0.03 for the ratio 
of the heavy-quark pole masses, we obtain: Xo+ (0) = 4.40^o;58 x 10"'^, ml Xi- (0) = 
9.92;[|:1? X 10-3, xo-(O) = 2.09l|j:[|^ X lO'^, and mfxi+iO) = 6.06+[J:i x lO'l For the 
spin-1 channels, we use the value mj, = (4.8 ± 0.2) GeV. The 0{as) corrections, which 



have been calculated in Refs. [^-[Q, enhance these results by 29%, 31%, 14%, and 
29%, respectively, where we use as = as{mb) = 0.22. The leading nonperturbative 



5 



corrections, proportional to the quark and gluon condensates, are very small and can be 
neglected [^, |2^. With these parameters, we obtain 

Xq+{0) = (5.69 ±0.71) X 10"^ 
mB-mD* Xi-(O) = (3.73 ± 0.51) x 10"^ , 

;^o_(0) = (1.46 ±0.06) X 10-2, 
mB-mD' Xi+ (0) = (3.63 ± 0.34) x 10"^ . (10) 

In the 1^ and 0^ channels, we have subtracted the pole contributions from the two 
lowest Be and B* mesons using the values of the decay constants shown in Table ^ 
In the evaluation of the inequalities (^, we choose to relate all form factors to the 
reference form factor Vi{w), which governs the B —>■ Dii? decay rate in (|I]), and whose 
normalization at zero recoil is V^i(l) = 0.98 ±0.07 Then, what enters the inequalities 
(|) are the products of the correlation functions divided by nfVi{l). The unitarity 
bounds become weaker the larger the values of the correlation functions and the smaller 
the value of Vi(l). Therefore, we simultaneously lower the value of the form factor 
and increase the central values in ([T0| ) by one standard deviation. This leads to the 
conservative upper bounds 

< 3.1 X 10-3 , ruB^mn^ < 2.0 x lO^^ , 



< 7.3x10-3, mB.mz5.^^^< 1.9x10-3, (11) 



n.Vni) ^ ' n,n2(l) 

which we use in our analysis. 



3 Model-independent bounds on form factors 

The four inequalities in (^ can be used to derive bounds on the corresponding form 
factors in the semileptonic domain. The problem is brought into a canonical form by 
performing a conformal mapping w z{w,a), which transforms the cut w plane onto 
the interior of the unit disc \z\ < 1, such that the integrals in (|^) become integrals around 
the unit circle z = e^^. This is achieved by defining 

z{w,a) = — a>0, (12) 



a 



which maps the branch point w = —1 onto z = —1. The choice of the real parameter a 
will be discussed later. 

As explained in previous works [T^-[T^, one further introduces a set of "outer" 
functions (f)j{z), i.e. functions without zeros nor singularities inside the unit disc, such 
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that the unitarity inequahties take the form 



271 

^JdeY.\Me'')F,H\" <l, (13) 

^ 

where it is imphed that w is expressed in terms of z = e*^ by means of (|T2]). Here Fj is one 
of the scalar, pseudoscalar, vector or axial- vector form factors, and j runs over the set of 
form factors in a given spin-parity channel. Along the unit circle, the modulus squared 
of the outer functions are equal to the weights in front of the corresponding form factors 
in (P), multiplied by the modulus of the Jacobian of the conformal transformation ([T^). 
The calculation of these functions is straightforward. The weight functions appearing in 
the inequalities are of the generic form 

A.(„^_l)3/.-..(!^)^"(,^_!^)-". (14) 

where iV is a constant. The corresponding outer functions are given by 

42-n 4 /I _i_^\3/2-n 




In deriving bounds on the values or derivatives of the form factors inside the unit disc, 
one must account for the effects of singularities located below the physical thresholds. 
For simplicity, we discuss the singularity structure in terms of the momentum transfer 
given in (^). For the form factors of interest, the thresholds occur at the values of \/q^ at 
which the relevant B^*^D^*^ pairs can begin to be produced, i.e. at {mB+mu) ~ 7.15 GeV, 
{niB + mD*) ^ 7.29 GeV, [mB* + nio) ~ 7.19 GeV, or (ms* + ttid*) ~ 7.33 GeV. Sub- 
threshold singularities arise due to the coupling to states with masses below these thresh- 
olds. On the one hand, there are cuts produced by multi-particle states. (Anomalous 
thresholds are seen not to appear by inspecting the relevant triangle diagrams at both 
quark and hadron level.) For instance, the form factor 5*3, which appears in the descrip- 
tion of B* D* transitions, has a branch point due to its coupling to BD pairs as well 



as to two-particle states consisting of a B^ meson and light hadrons [T^. However, these 
cuts are rather short. Their effects are suppressed by phase space and, in the case of Be 
states, by the Zweig rule. Once a model for the discontinuities across the cuts is adopted, 
these subthreshold singularities can be included in the dispersive bounds [0, and 
their effect is found to be very small. We thus neglect these contributions in the sequel. 

On the other hand, the form factors also have singularities due to single-particle 
states. Since both ground-state and orbitally excited Be states are expected to be very 
narrow ||2l|, we shall assume that they give poles located on the real axis. Our choice of 



form factors is such that the scalar functions Sj receive only contributions from scalar 
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particles, the vector functions Vj only from vector particles, etc., situated below the 
corresponding physical thresholds. The masses of these bounds states have already been 
given in Table |I[ The residues of the corresponding poles, which are proportional to 
the couplings of Be states to B^*^D^*^ pairs, are poorly known. Fortunately, dispersive 
bounds can still be derived even in the presence of poles with unknown residues 0. 
The optimal technique is to multiply each form factor Fj by a specific function Bj{z) (a 
so-called Blaschke or "inner" function) with zeros at the positions of the poles, and with 
unit modulus along the physical cut. The explicit form of the inner functions is 



n 



1 - ZZ^J^ 



(16) 



where the index n runs over all Be bound-states that couple to a given form factor, and 



(*) 



(17) 



is the image of the pole of mass M„ in the z plane. In terms of these functions, the 
unitarity inequality for a given spin-parity is written as 



< 1 



(18) 



where 



f,{z)=B,{z)4>,{z)F,{w{z)) 



(19) 



defines analytic functions inside the unit disc. We can thus Taylor expand the functions 
fj{z) in (0) and find that 



EE 

j n=0 



nl ■' 



< 1 



(20) 



where /j"^(0) denotes the n-th derivative of fj{z) with respect to z, evaluated at z = 0. 
There is one such inequality for each spin-parity channel. The derivatives with respect 
to z in ( PPD can be calculated using (|T^, where Bj{z) and (pjiz) are known functions. 
The derivatives of the form factors at 2; = are related to the derivatives with respect to 
w at the kinematical point w = Wq, defined such that z{wq, a) = 0. For the form factor 
Vi(w), for instance, we have 



dzViiwiz)) 
d'Mwiz)) 
d'Mwiz)) 



z=0 



z=0 



2 = 



-8aV?V^i(«^o), 

(128a% -32aV?) V'i(wo), 

(3072a^rfi + 1536a^ci - 144aV?) Vi{wo) 



(21) 
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where 



w{z) = 2a^ 



and the parameters pf, C\ and are defined by the expansion 

Vx{w) = Vi{wo) 1 - Pi(w - Wo) + ci{w - Wof + c?i(w - Wo)^ + 



(22) 



(23) 



Of course, the values of these parameters will depend on the choice of Wq, i.e. on the 
choice of the parameter a in (^). 

In the heavy-quark limit, all functions Fj{w) become identical in the semileptonic 
region and equal to the Isgur-Wise form factor. In order to incorporate corrections to 
that limit, we choose the vector function Vi{w) as a reference form factor and express 
the expansion parameters of all other form factors in terms of the reference parameters 
pI, c\ and d\ defined above. Our choice of the reference form factor is motivated by the 
fact that it is the physical form factor describing the semileptonic decay B ^ D iv . The 
leading symmetry-breaking corrections to the ratios of heavy-meson form factors have 
been analysed in great detail (for a review, see e.g. Ref. 0). We write 



w 



1 + Bj{w - Wo) + Cj{w - wo)^ + Dj{w - wof + ... 



(24) 



The results for the parameters y4j,...,Dj obtained by including the leading short- 
distance and l/mq corrections, for two different choices of a relevant to our discussion, 
are given in the Appendix. Using these results, the derivatives of the form factors Fj {w) 
can be expressed in terms of those of the reference form factor Vi{w). 



4 Zero-recoil expansion 

A convenient choice of the parameter a in (|I2p is to set a = 1 , in which case the conformal 
transformation w ^ z maps the zero-recoil point w = 1 onto the origin 2; = of the unit 
disc (i.e. wq = 1). With this choice, it is straightforward to obtain bounds on the form 
factors and their derivatives at zero recoil. We first give the allowed values for the slope 
pI and the curvature Ci, and then consider the coefficient di of the third-order term in 
(p3|). In order to see the effect of the corrections to the heavy-quark limit, we present the 
results both in the case of exact spin symmetry]^ and with the leading corrections taken 
into account. The effect of uncertainties in the corrections to the heavy-quark limit will 
be considered below. 

When only terms with up to two derivatives are kept in inequality (^0|) , the resulting 
constraint on the parameters pf and Ci of the reference form factor is given by an ellipse 
0: 

(P? - pD' + S Uc^ --c,)-T {pi - pD] < . (25) 



^In fact, we only equate the first three derivatives of the various form factors, but give all masses, 
poles and thresholds their physical values. 
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Table 2: Parameters of the ellipses for the different spin-parity channels, including 
corrections to the heavy-quark limit 





p1 
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K 




0+ 


0.58 


0.42 


1.08 


44.3 


1.05 


47.5° 


0- 


0.87 


0.81 


1.12 


61.9 


1.13 


48.4° 


1- 


0.86 


0.80 


1.14 


62.9 


3.84 


48.9° 


1+ 


0.95 


0.85 


1.42 


39.0 


3.92 


55.0° 



The single parameter sensitive to the values of the QCD correlators given in ( |Tl| ) is 
the "radius" K. All other parameters of the ellipse depend on known meson masses 
and, to a lesser extent, on the symmetry-breaking corrections of (^41). In Figure |l], we 



show the resulting ellipses for the different spin-parity channels, both with and without 
corrections to the heavy-quark limit. Only values of the parameters inside the ellipses 
are permitted. The effects of the corrections to the heavy-quark limit are small and 
tend to align the major axes of the four ellipses. In Table |^, we give the corresponding 
ellipse parameters, including the angle (p of the major axis. In all four cases, we find 
that S ~ 40-60, implying that the ellipses are highly degenerate. As a consequence, 
the curvature parameter ci is strongly correlated with the slope parameter pi, a model- 
independent property observed in Ref. |jT9|. It is apparent from Figure |I]that the bounds 



derived from the 0^ channels are much stronger than (and fully contained in) those de- 
rived from the 1^ channels. Therefore, we will not consider the latter hereafter. For 
completeness, we note that the (axial) vector and the (pseudo) scalar form factors are 
not fully independent functions, since some of them have to obey kinematical constraints 
at the point of maximum recoil, corresponding to = 0. When deriving the dispersive 
bounds in the different spin-parity channels, the kinematical constraints can be incorpo- 
rated using Lagrange multipliers. The dispersive bounds derived in the (pseudo) scalar 
channel are then interconnected with those derived in the (axial) vector channel, yielding 
a domain for the parameters of the reference form factor that is smaller than the one ob- 
tained from the intersection of the allowed domains obtained sperately from the different 
spin-parity channels. In practice, however, the (pseudo) scalar channel yields so much 
more restrictive bounds than the (axial) vector channel that the resulting improvement 
is negligible. 

To see what happens when the term with three derivatives is included, consider the 
simpler case of only a single form factor contributing to the unitarity sum. Then the 



corresponding inequality can be written in a form similar to (pSj): 

(pI - 6i)2 + (8ci - 62 - hplf + (64rfi - 64 - hpl - h^c^f < , (26) 
where the hi are numerical coefficients. This relation describes a highly degenerate 
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-4 -2 2 4 6 -4 -2 2 4 6 



2 2 
Pi Pi 

Figure 1: Bounds on the slope and curvature of the reference form factor Vi{w) 
obtained from the unitarity inequahties in the different spin-parity channels. Dashed 
curves correspond to the heavy-quark limit, while solid ones include the leading 
corrections to that limit. Note that the axes in the lower plots are extended by a 
factor of 2. 



ellipsoid in {pl,ci,di) space, whose axes are roughly in the ratios 1 : 1/8 : 1/64. The 
coefficient di is thus determined within a rather small range: 

d, = ^-^±M±^±A, (27) 
64 

where A < K/64 ^ 0.02. We will see below that this is much smaller than the range due 
to the theoretical uncertainty in the relations between the various form factors. What 
changes in the general case of many form factors are the precise values of the coefficients 
8 and 64 in front of Ci and di, as well as the values of the other ellipsoid parameters. 
In addition, the inclusion of the third derivative in ( pO] ) reduces the allowed domain 
for the ffist two derivatives. This feature is specific to the simultaneous treatment of 
several form factors whose derivatives are related, as they are here through heavy-quark 
symmetry. However, the inclusion of many form factors does not change the strong 
degeneracy of the eUipsoid in {pi, ci, di) space. 
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-0.5 0.5 1 1.5 2 2.5 -0.5 0.5 1 1.5 2 2.5 



2 2 
Pi Pi 

Figure 2: Bounds on the slope and curvature of the reference form factor Vi{w) 
obtained by including the third-derivative terms in the unitarity inequalities for the 
0^ channels, taking theoretical uncertainties in the calculation of corrections that 
break heavy-quark symmetry into account. 

Table 3: Parameters of the ellipses obtained from the third-order inequalities, in- 
cluding theoretical uncertainties. We also quote results for the third-order coefRcient 
di = api + (3ci + 7. 
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0+ 


0.68 


0.54 


1.01 


14.9 


0.85 


46.2° 


0.47(11) 


-1.33(6) 


0.01(11) 


0- 


0.86 


0.79 


1.11 


29.2 


1.08 


48.4° 


0.44(10) 


-1.37(6) 


0.05(11) 


Combined 


0.67 


0.55 


1.03 


25.9 


0.84 


46.4° 


0.45(7) 


-1.35(4) 


0.03(8) 



At first sight, the inclusion of the third derivative terms bring an improvement which 
appears to be very significant. For the case of the O"*" channel, for instance, the resulting 
allowed region for the parameters p\ and Ci is reduced, with respect to the one shown 
in Figure |l], by more than a factor of 3. However, the situation changes if we allow 
for theoretical errors 5Cj and 5Dj in the coefficients entering the relations (0) 
between the various form factors.^ To this end, we generate a large number of ellipses by 
taking random values of the theoretical errors within intervals [—5B, 5B\, [— 5C, 5(7], and 
[—5D, setting 5B = 6C = 6D = 0.1, which we believe is a conservative estimate of 
the uncertainties. The results are shown in Figure ^ Whereas the effect of theoretical 
uncertainties is moderate in the case of the second-order inequality, it has a large impact 
at third order. The envelope of the ellipses obtained in that case is much larger than each 
single ellipse, and almost fills up the ellipse obtained at second order. This shows that 

■^Tliere is no need to include errors in the normalization of the form fatcors at zero recoil, since we 
have already used a conservative low value for the normalization of the reference form factor. 
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1.5 

expansion in (w — 1) 
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w 




expansion in z{w,l) 



1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 
w 

Figure 3: Allowed shapes of the function for various values of p\. The dark 

(light) hands correspond to the third-order (second-order) expansion of the form 
factor. In the upper plot, the expansion is performed in powers of {w — 1), in the 
lower plot in powers of the variable z{w,l). 

not much can be gained by going to yet liigher orders in tlie expansion. Tlie tliick ellipses 
in Figure |^ are a good approximation to tlie allowed regions. The parameters of these 
ellipses are given in Table ^. Their intersection is again to excellent approximation given 
by an ellipse, whose parameters are given in the last row. We also present the results for 
the third-order expansion coefficient di, which we write as di = apl + (3ci + 'y. We find 
that the allowed range for this parameter obtained from (pTl) is only little affected by 
theoretical errors. The consistency between the results obtained from the two channels 
allows us to take the average of the two results, shown in the last row. 

To exhibit the possible behaviour of the form factor Vi{w) predicted by our results, 
we show in the upper plot in Figure ^ the allowed shapes for this function for a selection 
of equally spaced values of the slope parameter pi inside the allowed region, given by 
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—0.17 < p1 < 1.51. For each choice of pf, the parameters ci and di are scanned over 
the regions allowed by the dispersive bounds. The results are shown with and without 
including the third-order term in the zero-recoil expansion, i.e. the expansion in (^) 
with ti^o = 1. Close to zero recoil (i.e. for 1 < w < 1.2), both approximations provide an 
equally accurate description of the form factor. For larger recoil, the effect of the third- 
order term becomes clearly visible. The convergence rate of the series can be improved 
by expressing the form factor Vi{w{z)), with w{z) given in (^), as a power series in 
z rather than in {w — 1). The corresponding results are shown in the lower plot in 
Figure |^. Now the differences between the second- and third-order expansions become 
important only at larger values of w, and we are confident that the third-order results 
provide an adequate representation of the form factor over the enitre kinematical region 
accessible in semileptonic decays. Because the spread in the curves caused by variations 
of the parameter ci and di inside the allowed regions is very small, our results can be 
represented, to a very good approximation, by taking the values of these parameters 
along the major axes. This gives the one-parameter function 

Vi{w) 



1 - 8plz + {51.pl - I0.)z'^ - {252.pl - 84:.)z^ 



Vi{l) 

with . = ^^-f . (28) 

As an aside, we note that expanding this result in powers of {w — 1) would yield the 
slope-curvature relation ci ~ 1.05pl — 0.15, which is not very different from the relation 



ci ~ O.TApi — 0.09 obtained in Ref. [jT9[ with some assumptions about subthreshold 



singularities. We stress, however, that in the present work such assumptions are avoided. 



5 Optimized expansion 

We have seen above that higher-order terms in the expansion of the form factor Vi{w) 
around the zero- recoil point have a significant effect already for rather small values of w. 
The convergence of the expansion was improved considerably by introducing the variable 
z instead of {w — 1); however, the question arises whether an even better convergence 
can be achieved by optimizing the choice of the expansion variable. 

When one wishes to approximate a function along an interval, the most natural 
expansion is given in terms of orthonormal polynomials on this interval. The domain 
of convergence of the expansion is then an ellipse passing through the first singularity 
of the function ||2^. In order to accelerate the rate of convergence in the interval of 



interest, one must conformally map the interior of the function's analyticity domain 
onto the inside of an ellipse, such that the interval of interest is applied to the segment 



between the focal points [^. This conformal mapping is given by elliptic functions 
P^ . However, when the physical interval is far from the first singularity (as it is for 
semileptonic -B*-*^ D^*^ decays), the ellipse is almost degenerate to a circle, so that 
a simple Taylor expansion about the center of the ellipse will have almost the same 
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radius of convergence as the expansion in orthonormal polynomials discussed above. In 
the absence of subthreshold singularities, this approximate circle is simply a rescaling of 



the circle obtained using the conformal mapping of (12), with a determined such that 
the physical region is applied onto an interval symmetric about the origin z = 0. This 
gives a = Qopt such that 2;(wmax, ^opt) = —z{l,aopt), which is very close to the variable 



first suggested by Boyd and Lebed [^. For the case oi B ^ D transitions, we find 



ttopt ^ 1.067 and z^ax = ^(wmax, ctopt) ~ 0.032. 

The argument just presented is only valid when there are no subthreshold singular- 
ities, while the form factors of interest have poles and branch points below the physi- 
cal cut. Nevertheless, because the singularities are far from the semileptonic domain, 
the variable z{w, Copt) remains an excellent choice. Indeed, a measure of the rate of 
convergence of an expansion of a form factor in a variable x is provided by the ratio 
'^conv = |a;|max/|a^poie| whcrc x^oXe IS the position of the first pole and |x|max the largest 
value can take in the semileptonic domain. The form factor Si{w), for instance, has 
a pole due to a scalar state at 6.70 GeV, so that rconv ~ 0.07. For comparison, we note 
that Tconv ~ 0.35 for the expansion variable x = {w — 1), and rconv ~ 0.15 for the variable 
X = z{w, 1) used in the previous section. It is clear that z{w, Oopt) gives the best rate of 
convergence. In order to simplify notation in what follows, we will write z^^ for z{w, Oopt) 
and a* for Oopt- We shall use the same value of a* for all form factors. 

Our claim is that all B^*^ D^*^ form factors can be accurately described by a 
second-order polynomial in z^,. In particular, we shall derive dispersive constraints on 
the slope and the curvature Cu (now taken at the recoil point w = Wq ^ 1.276) in 
the expansion of the reference form factor Vi{w{z^:)), i.e. 

^'^"^^ - 1 - Salplz, + 16al{Aalcu - pl)zl . (29) 



Viiwo) 



To estimate the truncation error due to the neglect of higher-order terms, consider the 
expansion of the analytic functions fj{z) of (|19D around z = 0. If we write 

M.)^MO){i^^^^^j^M^)). (30) 

the remainder Aj{z) can be bounded using the dispersive constraint (^0]). We find that 



l^-WI<l^iy^^^-l- (31) 



Let us concentrate, for concreteness, on the scalar form factor Si{w). Making the con- 
servative assumption that Si{wo) is as low as |S'i(l), we find |A5^(2;)| < 7 x 10"'^ for 
all values of z corresponding to the semileptonic region. This truncation error is tiny 
and will be neglected it in what follows. The second step consists in noticing that the 
known function [(f)Si{z^:)Bs-i^{z^)]~^ is well approximated by a second-order polynomial 
in z*. This is where the convergence criterion discussed earlier comes into play, since 
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Figure 4: Bounds on the slope and curvature of the reference form factor Vi{w) 
obtained using an expansion in the optimized variable z*. The solid (dashed) ellipses 
show the results obtained with (without) including theoretical uncertainties in the 
heavy-quark relations. 



l/i?5^(-2*) contains the scalar poles. The fact that rconv = k*|max/k*poie| is so small 
guarantees that the expansion of [05^(z,,)i?5^(z,,)]~^ converges rapidly. We find that the 
first three terms in this expansion deviate by at most 0.2% from the full function over 
the whole semileptonic region. As a result, the form factor Si{w), which is proportional 
to the product of the two quantities fsA^^) and [(t)s-i{.z^)Bsi(z^)]~^ , is approximated to 
an accuracy of better than about 0.2% by the product of the truncated expansions of 
these two factors. Of course, this product includes terms of order and z'^. However, 
their contribution never exceeds 1.5% over the full semileptonic region (it is lower than 
1% in the region w < 1.5). Thus, we conclude that the scalar form factor Si{w) can be 
described by a second-order polynomial in z^, with an accuracy of better than 2%. A 
similar conclusion holds for all other B D'^*^ form factors, since they are related to 
the Si{w) by heavy-quark symmetry, up to small corrections. 

As in Section ^, we use the dispersive bounds to constrain the first two derivatives 
of Vi{w) at wq. However, wq is not the zero-recoil point here, and the left-hand side of 
(PPP contains as an overall factor the unknown normalization |Vi(wo)P- To proceed, we 
multiply both sides of the inequality by a common factor. 



EE 

j n=0 



< 



Viiwo) 



(32) 



such that now the left-hand side contains the known normalization Vi(l). To evaluate 
the right-hand side of (|3^ ) we rely on our observation that in the semileptonic region 
all D^*^ form factors are very well described by a second-order polynomial in 

z*, and thus use (p9[) with w = 1. This approximation enables us to write the resulting 
constraints on p^^ and Ci^, as ellipses, like in the zero-recoil case. 

For the scalar and pseudoscalar channels, the ellipses obtained in the (p^*, ci*) plane 
are shown in Figure The dashed ones correspond to using the central values of the 
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Table 4: Parameters of the ellipses obtained from the second-order inequalities, 
using the optimized variable 





-pl 


Cl 


T 


S 


K 




0+ 


1.11 


0.89 


0.90 


29.0 


1.40 


42.5° 


0- 


1.07 


0.88 


0.89 


27.6 


1.26 


42.2° 



parameters Aj, Bj and Cj defined in and given in the Appendix. As in Section ^ 
we account for tlieoretical uncertainties in tlie values of these parameters by generating 
a large number of ellipses, allowing Bj and Cj to vary by ±0.1. This procedure generates 
the shaded areas, which to a good approximation are contained in the larger ellipses. We 
find again a similarity between the final ellipses in the two channels and, in particular, a 
near alignement of their major axes. The intersection of these two ellipses is, to a very 
good approximation, given by the pseudoscalar ellipse alone. This is in contrast to the 
zero-recoil expansion considered in the previous section, where it was the scalar channel 
that gave the tighter constraints. For completeness, the ellipse parameters are given in 
Table |. 

Let us now investigate the properties of the optimized parametrization and compare 
them with the results obtained in the previous section. In the upper plot in Figure ^ 
we show the allowed shapes of the form factor Vi{w) for values of the slope parameter 
pI^, chosen such that they correspond to the values of the zero-recoil slope pl used in 
Figure For each choice of p^^, the parameter cu is scanned over the region allowed 
by the dispersive bounds. Because S is large compared with K^, and because the term 
multiplying cu in expression (^9]) for the form factor is suppressed by the small factor z^, 
the spread in the curves caused by the variation of cu is very small. Therefore, to a very 
good accuracy, Cu can be replaced by ci + tan ip {pl^ — p\). This gives the one-parameter 
function 

^ l-9.105pi^. + (57.0pL-7.5)z,^ 
Vi{l) ^ 0.354pi + 0.992 

with = ^ . . (33) 

v/wTT+ 1.509 ^ ' 

The light-shaded bands in the figure include, in addition to the spread discussed above, 
an estimate of the truncation error due to the neglect of higher-order terms in z*. It 
is obtained by not expanding the known functions \^j{z^)Bj{z^~^ in powers of z^, 
but keeping their full expressions to calculate the form factors from the second-order 
expansions of the functions (We have argued above that the error from the 

truncation of the functions /j(-2*) themselves is a negligible effect.) These bands reflect 
the total theoretical uncertainty in our results, which is seen never to exceed the level 
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Figure 5: Upper plot: Allowed shapes for the function Vi{w) for various values of 
Pi^ (the slope at wq « 1.276) chosen such that they correspond to the values of 
Pi (the slope at w = 1) of Figure |^. The light bands include an estimate of the 
truncation error. Lower plot: Comparison of the two one-parameter approximations 
(H) (solid) and @ (dashed). 

of 2%. As a consequence, the one-parameter description of the form factor given above 
can be trusted to that level of precision. 

In the lower plot in Figure ^, we compare the parametrization ( PBf ) to the one given 
in (^). The parameters pi and are chosen as before. We find that the two sets of 
curves are nearly indistinguishable in the semileptonic region 1 < w < 1.6, whereas they 
start to diverge for larger w values. This agreement provides a strong consistency check 
of our method, given that the expansions in the two variables z and 2* are carried out 
to different order, and that the dispersive bounds in the two cases are dominated by 
different spin-parity channels. We thus conclude that, for all practical purposes, both 
parametrizations are equally reliable. Similar parametrizations for all other B'^*^ D^*^ 
form factors can be obtained using p^ ) and the results given in the Appendix. 
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6 Phenomenological applications 



We now summarize our results and provide simple, constrained parametrizations for the 
form factors which determine the differential rates for B ^ D iu and B — > D*i v decays. 
The form factor Vi [w] governing the first process is described to an accuracy of better 



than 2% by the two parametrizations given in (|2l 
entering these expressions are related by 



and (^3|) . The two slope parameters 



Pi* 



0.62p2 + 0.04 



1.61p?, - 0.06 
1 + 0.36pL 



(34) 



1 - 0.22p? ' 

The dispersive bounds discussed in Section || require that —0.17 < pi < 1.51, and 
therefore —0.07 < pf^ < 1.47. These intervals, by themselves, represent nontrivial 
bounds on the slope of the form factor, which are competitive with the most recent 
bounds derived from inclusive heavy-quark sum rules [ pij] . 

To obtain similar parametrizations of the function T{w), which governs the B — > 
D*iu decay rate in (0), we find it convenient to first relate this function to the axial- 
vector form factor Ai{w) through 



1 + 



Aw 1 — 2wr -\- r 



'1-r 



1 — 2wr + r 



:i-r) 



w 



1 



w + l 



+ 



1 + 



w 



Ai{wf, (35) 



where r = mo'/^Bi and Ri{w) and R2{w) are ratios of form factors defined in the 
Appendix, and are given by 

hv{w) 



Riiw) 



R2(W) 



w 



1.27 - 0.12{w - 1) + 0.05(w; - 1) 



» + rhA^w) 



0.80 + 0.11(w - 1) - OM{w - 1) 



(36) 



The theoretical predictions for these ratios are supported by measurements reported by 
the CLEO Collaboration ||32|. We have chosen to work with Ai{w) instead of J^{w), 
because the latter suffers from large kinematical enhancements of the corrections to 
the heavy-quark limit. In other words, the coefficients in the expansion of the ratio 
J-'{w)/Vi{w) corresponding to (|37D below would be large. These large corrections were 
not accounted for in Ref. where the symmetry-breaking corrections in the relation 
between J-'{w) and Ai{w) are neglected. 

Using the results of the Appendix, the form factor Ai{w) can be related to the 
reference form factor Vi{w) through 



Vi(w) 



0.948(l - 0.212Z - 4.007^2 - 1.342z^ + 
0.476^,-4.163^2 + ...' 



0.937(1 



(37) 
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Table 5: Results of the various fits to experimental data on the recoil spectra in 
B — > D*ir' decays j |3^ . The slope parameters are measured at w = 1. 



Parametrization 








Eq. ® 


0.037(2) 


1.36(18) 


1.15(18) 


Eq. (H 


0.037(2) 


1.34(17) 


1.13(17) 


Linear Fit 


0.035(2) 




0.81(13) 



where z = z{w,l) and z^, = z{w,a^:). From this, it is straightforward to derive the 
parametrizations of Ai{w) from the results for Vi{w) given in (|28|) and (0). They are 

^ 1 - Spl^z + {53.p\, - 15.)^' - (231.pi, - 91.)z' , 
with . = ^^-f , (38) 

and 

A,{w) _ 1 - 9.105pX,z, + (61.3pi^, - 14.8)z,^ 
Ai(l) 0.358p^^, + 0.984 

with z^: = . (39) 

VwTl + 1.509 ^ ' 

The relation between the two slope parameters entering these expressions is 

^ O.eOpj^ + 0.07 l-66pi,.-0.12 

l-0.22pi^ ' ^"-'^ l + 0.36p^^, ■ ^ ^ 

The dispersive bounds discussed in Section ^ require that —0.14 < < 1.54 and 
—0.01 < p\_^^ < 1.51. Alternatively, we could repeat the derivation of the dispersive 
bounds using Ai{w) instead of Vi{w) as the reference form factor. Since the symmetry- 
breaking corrections in the relation between the two form factors are very small, the two 
procedures give essentially the same results. 

Our results will eliminate the large uncertainty associated with the extrapolation of 
the experimental recoil spectra in i? ^ D^*^ii? decays to zero recoil. At present, this 
uncertainty is the main experimental error in the determination of |Vc6| [Q. In fact, 
in many experimental analyses a linear shape of the form factor has been used to ex- 
trapolate the data, without attributing an error to this assumption. A reanalysis using 
our constrained parametrizations is therefore highly desirable in order to get a reliable 
determination of \Vcb\- To illustrate the potential impact of our results, we show in 
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Figure 6: Experimental data for the product {Vcbl^iw) as a function of w, ex- 
tracted from semileptonic B D*£iy decays J3^/ . The sohd curve is a fit using 
our parametrization in while the (barely visible) long-dashed one refers to the 



parametrization of (3£ ). The two curves are nearly indistinguishable in the semilep- 
tonic domain. The short-dashed line shows a straight-line fit, similar to the one used 
in the experimental analysis. 



Figure ^ different fits to the data for the product J-'{w) \ Vcb\ measured in B ^ D* 



decays by the CLEO Collaboration |Q. Similar measurements have also been reported 
in Refs. |^-|3^. Our goal here is not to establish a new value of \Vcb\ (which is best 
done by the experimental collaborations themselves), but rather to illustrate how our 
parametrizations work. The short-dashed line shows a linear fit to the data, whereas the 
solid and long-dashed curves show the fits obtained using our parametrizations. All fits 
have an excellent x^/dof between 0.3 and 0.4, reflecting that we are neglecting corrre- 
lations between the errors of the various data points. Clearly, the two parametrizations 
we propose give nearly indistinguishable fits. More importantly, the obtained values of 
\Vcb\ and the zero-recoil slope parameters are systematically larger than those derived 
from the linear fits. 

The results of the fits are summarized in Table ^, where we also quote values for the 
slopes of the form factors Ai{w) and J^{w) at zero recoil. The theoretical prediction for 
the difference between the two slope parameters, which can be derived using the results 
of the Appendix, is p% ^ - 0.21 0. 



7 Conclusions 

We have derived dispersive bounds on the shape of the form factors describing semilep- 
tonic B D^*Hiy decays. The method we use is based on dispersion relations and 
complex-analysis techniques, which combine QCD calculations of current-current corre- 
lation functions in the Euclidean domain with spectral representations of these functions 
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in terms of sums over intermediate hadronic states. Using the positivity of these spectral 
functions and crossing symmetry, it is possible to derive bounds on the form factors of 
particular hadron states that are contained in the spectral sum. In the case of heavy 
mesons, these bounds can be strengthened by including in the unitarity sum all contribu- 
tions that are related, in the semileptonic region, by heavy-quark spin symmetry. Thus, 
we include the contributions of BD, BD*, B*D and B*D* states, taking into account 
the leading short- distance and l/mg corrections to the heavy-quark limit. We consider 
all four spin-parity channels relevant for B^*^ D*^*^ transitions and combine them in 
an optimal way. 

We first apply these methods to derive bounds on the slope, curvature and third 
derivative of a reference form factor in an expansion around the zero- recoil point w = 1, 
using as input the value of this form factor at zero recoil. We find that a consistent 
inclusion of the third-derivative term requires a careful treatment of theoretical uncer- 
tainties in the relations that connect various form factors in the heavy-quark limit. We 
show explicitely how accounting for these uncertainties reduces significantly the apparent 
improvement brought by the third-derivative term. 

In analysing the different spin-parity channels, we find that the scalar and pseu- 
doscalar channels give significantly better bounds than the vector and axial-vector ones. 
We find a strong correlation between the curvature and slope of the reference form factor, 
and an even stronger correlation of the third derivative with the slope and curvature. 
This observation has led us to investigate expansions in other kinematical variables, in 
which the strong correlations between the parameters of the {w — 1) expansion are ex- 
plained naturally. We find that the variable z{w,a) defined in ( ]T^ ) accomplishes this 
task, and furthermore leads to an improved convergence of the expansion of the form fac- 
tors. We have discussed two choices of the parameter a that are particularly convenient: 
the choice a = 1 leads to an expansion around zero recoil, however in a better variable 
than (w — l); the choice a = discussed in Section || corresponds to an expansion around 
a kinematical point close to the center of the semileptonic region, thus providing nearly 
optimal convergence. For the first choice, we find that the scalar channel provides the 
tightest constraints, and that the form factors can be well approximated by third-order 
polynomials. For the second choice, the pseudoscalar channel gives the best bounds, and 
the form factors can be approximated by second-order polynomials. In both cases, we 
derive simple one-parameter functions for the physical form factors Vi{w) and Ai{w), 
which are accurate to better than 2% in the semileptonic region. Similar parametriza- 
tions for all other 5^*^ D^*^ form factors can be obtained using (p^) and the results 
given in the Appendix. The fact that the parametrizations obtained in the two cases 
lead to shapes of the form factors that are nearly indistinguishable in the semileptonic 
region provides a strong consistency check of the method. We note that our functions 
are simpler than those derived by Boyd et al. [^] and include the leading short-distance 
and l/friQ corrections to the heavy-quark limit, as well as the uncertainties in their cal- 
culation. Nevertheless, we have checked that in fits to available experimental data our 
results for the form factors agree with those obtained in Ref. pO[ to better than 3% 
for values w < 1.5. For larger values of w, the parametrization for T{w) given in this 
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reference becomes less reliable. 

For the convenience of the reader, we restate our results obtained in the zero-recoil 
expansion based on the variable z{w, 1), which gives the simplest parametrizations. They 
are 



Vi{w) 





(1) 




w) 




(1) 



1 - Spfz + {51.pi - 10.)^ - (252.pf - 84.)^^ , 
1 - Spl^z + {53.pl^ - 15.)z' - (231.p^^ - 91.)^' , (41) 

where z = (^/w + 1 — \/2)/{\/w + 1 + -\/2), and p\ and p\_^ are the slope parameters 
at zero recoil, restricted to the intervals —0.17 < p\ < 1.51 and —0.14 < p\^ < 1.54. 
The values of the form factors at zero recoil are known to be Vi(l) = 0.98 ± 0.07 
and Ai{l) = 0.91 ± 0.03. In Section |^, we have shown how these results can be used to 
determine \ Vcb\ from an extrapolation of experimental data on semileptonic decays to zero 
recoil. Not only do our parametrizations eliminate the main experimental uncertainty 
in this extraction, but they also lead to a corrected value for |Vcfc| that is systematically 
higher (by about one standard deviation) than the value obtained from a linear fit. 
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Appendix: Form factors and heavy-quark relations 



In the context of the heavy-quark expansion, it is common practice to parametrize the 
meson matrix elements of the weak currents = c'y'^b and = C7'^7^6 by a set of 
form factors hi{w) that depend on the kinematical variable w = v ■ v', where v and v' 
are the meson velocities. With a mass-independent normalization of meson states, these 
functions are defined as 

{D{v')\V^'\B{v)) = h+{w) [v + vY + h^w) [v - v'Y , 

{D*{v',e')\V''\B{v)) = 2/iv(«;)e''"°^£»/3, 

{Div')\V''\B*iv,e)) = ihy{w)t^^''^e,v'^vp, 



{D*{v\e')\V^\B*{v,e)) 



hi{w) {v + v'Y + h2{w) {v - v'Y 
+ h{w) e'* ■ V e^" + h^iw) e ■ v' e'*" 



/* 

e ■ e 



I* I 
e ■ ve ■ V . 



(A.l) 



for the vector current, and 

{D\v\e')\A^\B{v)) 

{D{v')\A^\B\v,e)) 
{D*{v\e')\A^\B*{v,e)) 



hA, (w) {w + 1) e'*^ - [hA, {w) v^" + hA, (w) v'f" 
hA, (w) {w + 1)6^"- \hA, (w) v''' + Ha, (w) v'' 



e ■ V . 



e ■ V 



le 



+ 



hg{w) e'* ■ vSy + hio{w) e ■ v' e 



I ^1* 

V 



(A.2) 



for the axial-vector current. Here, e'-'-* denote the polarization of the vector mesons. 
For our purposes, it is convenient to introduce a new set of form factors defined as 
linear combinations of the functions hi{w). We introduce these form factors in such a 
way that they reduce to the Isgur-Wise function in the heavy-quark limit, have definite 
spin-parity quantum numbers, and "diagonalize" the the unitarity relations derived in 
Section |. Thus, we define the scalar ( = 0"*") functions (dropping the argument w for 
brevity) 



(BD) 



g{B*D*) 



1 + r w — 1 

1 — r w + 1 



1 + r w — 1 

hi - — h2 , 

I — r w + I 

1 + r w — 1 



w ( hi 



+ 



1 — r w + 1 

W — 1 r 



1 — r 



r — hi + {1 — wr) h5+ {w — r) ht 



(A.3) 
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the pseudoscalar (J^ — ) functions 



PI 



(BD") 



y{B*D*) 



1 + r 

1 

1 + r 
h7- 



[w + 1) /iAi - (1 - wr) — {w - r) Haz 



r{w + 1) hA, -{r-w) Ha^ - {rw - 1) Ha^ 



1-r 



1 + r 



the vector (J^ — 1 ) functions 



(ED) 



h+ — h- , 



{B*D*) 



1+r 
1 ~ r 
1+r 



h2 , 



' ^w(hi- — — /is 
1 + r 





+ 




= hv 


y(B*D) 


= hv 






yiB*Dn 


= hi 



1 + r 



(1 - wr) h3 + {r - w) h^ + (w^ - 1) {rh^ + /ie) 



and the axial- vector (J^ = 1+) functions 



^3 



hAi , 

W — 1 

h'j - , , /is + (ui - 1) /iio , 



/i7 + 

1 

1 — r 
1 

1-r 



w + 1 
w — 1 
+ 1 



/is + - 1) /ig , 



(w - r) /iAi - - 1) (r/iA2 + /i^s) 
(1 - wr) hA, + (w-l) (hA^ + rhAs) 



AB*D*) , 1 + rw-l 

A\ ' ^ hr- — - /is 

1 — r + 1 
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In these equations, r = mj^(,) /m^i,) denotes the appropriate ratio of meson masses. 
The functions hy and h^., which have not been defined in Ref. 0, are obtained from 
the corresponding functions hy and /ia^ by an interchange of the heavy-quark masses, 
rric ruh- 

Luke's theorem protects the functions /;.+, /iai, ^Ai? ^i? ^ind hj from receiving first- 
order power corrections at zero recoil Up to corrections of order l/mg, they obey 
the normahzation conditions 



S,{l)=7^V, A.j{l) = 7^A, (A.7) 

where rjv and t^a are the QCD renormahzation constants of the vector and axial cur- 
rents at zero recoil [^. This implies that, to order l/rriQ, the functions appearing in 
the unitarity inequalities for the spin-parity channels O"*" and l"*" in @ have the same 
normalization. 

In the heavy-quark limit, h+ = hy = /^Ai = = hy = h^^ = h^.^ = hi = = 
hi = hj = ^, where ^{w) is the universal Isgur-Wise function, while all other functions 
hi vanish. The leading short-distance corrections to this limit can be parametrized in 
terms of the Wilson coefficients Cj and Cf of the vector and axial currents {i = 1, 2, 3), 
which are calculable functions of the variable w and the heavy-quark masses, and 
rric- For our purposes, it is sufficient to work with the exact one-loop expressions for 
these functions [^-[^. Leading logarithms, which have been summed to all orders in 



perturbation theory, are universal and cancel in the form-factor ratios we are interested 
in. In fact, for the same reason we are free to divide out Ci as a common factor. Then 
the resulting one-loop expressions are 



Cf^ 2«, 



H(5)iw, 1/z) , 



T^H^5){w,z), (A. 



Ci 3n 
where z = nic/mb, and 

r{w) = . In (w + Vw'^ — l) , 

— 1 

^(5)(^, = f~^'"'r2 + n o\ 2,2 |2(^ T mi ± z) In. 
l-2wz + z'^ {I -2wz + z^y I 

- \w ± 1) - 2W{2W ± 1)Z + (5W ± 2^2 ^ 1)^2 _ 2^3j ^(^y^^^j ^ (^^g^. 

The upper signs refer to the function the lower ones to H^. 
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Similarly, the leading power corrections can be parametrized by six functions Li{w) 
introduced in Ref. |^, which are related to the subleading Isgur-Wise functions ri{w) and 
Xi(u') defined in Ref. 0]. Once again, since we are only interested in form-factor ratios, 
we can omit the contribution of xii'w)^ which is the same for all meson form factors. 
Then then resulting expressions are 

Li = -4{w - 1) ^ + 12 ^ ^ 0.72(w - 1)A , 
L2 = -4 y ^ -0.16(w - 1)A , 

L3 = 4 ^ ^ -0.24A , 

L4 = (2r7 - 1) A ^ 0.24A, 
L5 = -A , 

2 - .S 94 _ 

L, = — (ry + 1) A ^ A , (A.IO) 

w + 1 w + I 

where A = rriM — ttlq 0.5 GeV is the "binding energy" of a heavy meson, and on 
the right-hand side we have used approximate expressions for the subleading Isgur-Wise 
functions obtained using QCD sum rules [Q. 

The corresponding expressions for the functions hi are 

y = C7i + (C2 + C3) + {Sc + Sb) Li , 

y = (C2 - C3) + (e, - Et) L4 , 

hy 



^ = Ci + Ec — , ^ L5) + Eb (Li — , ^ I , 

^^2 _ ^5 



A3 



T 



Cf + C| + {L2 -L3-L5 + Le) + (Li - L4) 



Ci + {Li - L4) + Eb {L2 - L 



5J , 



^Ai _ ^5 W - 1 , \ W-l 



CI + Ec [Li — LA+Eb[L2 — L5 



^ \ w + 1 J \ w 
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i 
h 



-CI + £b (L3 + Le) , 
-^ = Cl- Cl + (Li - L4) + (L2 - L3 - L5 + Le) 
^ = Ci + ^(C2 + C3) + (ee + e,)L2, 

/l4 

T 



Cl + [L2 + - 1) L3 + L5 - (w + 1) Le] + Sb {L2 - L^) 
Cl + e, (L2 - L5) + [L2 + - 1) L3 + L5 - + 1) Le] 

-C2 + £c (^3 - L^) , 
— C3 + £b {L3 — Lq) , 



hs w + l 5 



{C',+Cl) + {e,-eb)L,, 



e 2 

C2^+£, (L3-L6) 



Cl + eb{L,-U), (A.ll) 



where = l/2mb and = l/2mc. Given these exphcit expressions, it is straightforward 
to evaluate the form-factor ratios defined in p^ ) including the leading corrections to the 
heavy-quark limit. In order to be consistent, we expand the expressions for the ratios 
to linear order in and l/mq. In Table [A.l| , we show the results for the expansion 
coefficients in (^If), with wq = 1. In the numerical analysis, we use = as{^/mijn^) = 
0.26, z = mc/rrih = 0.29, and A = 0.48 GeV. In Table |A.2| , we give the corresponding 
values for the same expansion parameters, but at Wq ~ 1.267. 
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Table A.l: Theoretical predictions for the coefRcients in the expansion of the form- 
factor ratios Rj (w) around wq = 1 





A, 


B, 


C, 




Si 


1.0036 


-0.0068 


0.0017 


-0.0013 


S2 


1.0036 


-0.0355 


-0.0813 


0.0402 


S3 


1.0036 


0.0776 


-0.1644 


0.0817 


Pi 


1.1548 


-0.2088 


0.0032 


-0.0009 


P2 


0.9060 


-0.0727 


0.0031 


-0.0008 


P3 


1.0325 


-0.2116 


0.0032 


-0.0009 


Vi 


1 











V2 


1.0681 


-0.1944 


0.0000 


0.0000 


V3 


1.1361 


-0.2474 


0.0000 


0.0000 




1.2179 


-0.1428 


-0.0015 


0.0004 




1.0676 


-0.0362 


-0.0015 


0.0004 




1.4919 


-0.2278 


-0.0015 


0.0004 


Vi 


1.3416 


-0.1987 


-0.0015 


0.0004 


Ai 


0.9484 


-0.0265 


-0.0560 


0.0266 


A2 


0.9484 


0.0050 


-0.0184 


0.0078 


A3 


0.9484 


-0.0205 


-0.0869 


0.0421 




0.9484 


0.0256 


-0.1245 


0.0609 


A, 


0.9484 


0.2984 


-0.2391 


0.1180 


^6 


0.9484 


-0.1699 


0.0621 


-0.0324 


A, 


0.9484 


-0.0113 


-0.0857 


0.0414 
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Table A. 2: Theoretical predictions for the coefEcients in the expansion of some of 
the form-factor ratios Rj{w) around wq ^ 1.267 



F- 

^3 


A ■ 






Si 


1.0018 


-0.0061 


0.0009 




0.9883 


-0.0727 


-0.0554 


Ss 


1.0140 


0.0025 


-0.1117 


Pi 


1.0974 


-0.2072 


0.0026 


P2 


0.8862 


-0.0711 


0.0025 


Ps 


0.9743 


-0.2100 


0.0025 


Ai 


0.9373 


-0.0522 


-0.0387 
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